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Abstract 

We study the chiral phase transition at finite temperature T and baryonic chem- 
ical potential within the framework of the SU(3) Nambu- Jona-Lasinio (NJL) 
model. The QCD critical end point (CEP) and the critical line at finite T and 
[ib are investigated: the study of physical quantities, such as the baryon number 
susceptibility and the specific heat in the vicinity the CEP, will provide relevant 
information concerning the order of the phase transition. The class of the CEP is 
determined by calculating the critical exponents of those quantities. 
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1 Introduction 



It is commonly accepted that the vacuum of quantum chromodynamics (QCD) 
undergoes a phase transition to the quark gluon plasma (QGP) at high tem- 
perature and/or quark chemical potential. Such a new state of matter is ex- 
perimentally studied in on-going heavy ion collisions at CERN, Brookhaven 
and JIRN pQ. 
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The discussion about the existence of a tricritical point (TCP) or a CEP is 
also a topic of recent interest. As is well known, a TCP separates the first 
order transition at high chemical potential from the second order transition 
at high temperatures. If the second order transition is replaced by a smooth 
crossover, a CEP which separates the two lines is found. The existence of the 
CEP in QCD was suggested at the end of the eighties [2|3] . and its properties 
have been studied since then (for a review see Ref. [HE]). 

The most recent lattice results with Nf = 2 + 1 staggered quarks of physical 
masses indicate the location of the CEP at T CEP = 162±2MeV, yF EP = 360± 
40 MeV |6] , however its exact location is not yet known (it depends strongly of 
the mass of the strange quark). At the CEP the phase transition is of second 
order, belonging to the three-dimensional Ising universality class, and this 
kind of phase transitions are characterized by long-wavelength fluctuations of 
the order parameter. 

The possible signatures of the CEP in heavy-ion collisions have been studied 
in detail in [4"]TT]|8] . In heavy-ion collision experiments, fluctuations of T and 
fiB, can be found, respectively, in event-by-event fluctuations of pt spectra or 
in event-by-event fluctuations in the baryon number to pion ratio [Tll§] . The 
event-by-event fluctuations of T can be related to the heat capacity: if the 
specific heat C diverges, as the CEP is approached from either the left or the 
right, the fluctuations of T are suppressed. On the other hand, the event-by- 
event fluctuations of \lb can be related to the baryon number susceptibility, xb, 
which is also divergent. This implies that fluctuations of \iq are also suppressed 
at the critical pointf"*"!. 

As pointed out in Ref. [10] , the critical region around the CEP is not pointlike 
but has a richer structure. The critical region is defined as the region where 
the mean field theory of phase transitions breaks down and nontrivial criti- 
cal exponents emerge. The size of this critical region is important for future 
searches for the CEP in heavy ion-collisions |llj . 

It is also expected that the strange quark should have an important effect on 
the position of the TCP and the CEP. At finite T and zero /ig, in the limiting 
case with m u = = and infinite strange quark mass m s , the chiral phase 
transition is likely to be of second order and the static critical behaviour is 
expected to belong to the universality class of the Heisenberg 0(4) model in 
three dimensions [12]. When m s is finite and less than some critical value mf lt , 
the second order transition becomes of first order. This leads to a tricritical 



The divergence of xb is directly related to an anomaly in the event-by-event 
fluctuation of baryon number: in a heavy-ion collision experiment it is expected 
that the event-by-event fluctuation of the proton number is relatively enhanced 
for collisions which have passed in the vicinity of the CEP or the TCP. It is also 
expected an increase in event-by-event fluctuations at low pt near the CEP |4|7j . 
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point in the T — m s plane [13| . 



Some studies have been done in the SU(2) sector [TUfTi] but less attention has 
been given to the effects of the strange quark [15]. In this paper we aim at 
investigating the chiral phase transition and the CEP in quark matter with 
strange quarks. 



2 Model calculations 



We perform our calculations in the framework of the three-flavor NJL model 
prB1H7pi8lTTU] , including the determinantal 't Hooft interaction that breaks the 
£7^4 (1) symmetry, which Lagrangian reads 



£ = q (id ■ 7 - m) Q + f £ [(<?A"g) 2 + (q(i l5 ) \ a qf 



^ a=0~ 



+ 9D 



det[g(l + 75 )g] +det[g(l- 75 )g] 



Here q = (u, d, s) is the quark field with three flavors, Nf = 3, and three colors, 
N c = 3. rh = diag(m u , m^, m s ) is the current quark mass matrix and A a are the 
Gell-Mann matrices, a = 0, 1, . . . , 8, A = I. The three momentum integrals 
are regularized by the cutoff A, and a standard set of parameters [TTf2"0] given 
by A = 602.3 MeV, g s A 2 = 3.67, g D A 5 = -12.36, m u = m d = 5.5 MeV 
and m s = 140.7 MeV allow to reproduce the following vacuum observables: 
M w = 135.0 MeV, M K = 497.7 MeV, / ff = 92.4 MeV and M v > = 960.8 MeV. 
With this set of parameters we also obtain (q u q u ) = (qdqd) = —(241.9 MeV) 3 
and (q s q s ) = — (257.7MeV) 3 , for the quark condensates and M u = Md = 367.7 
MeV and M s = 549.5 MeV, for the constituent quark masses. 

The fundamental relation is provided by the baryonic thermodynamic poten- 
tial 

n(jH,T) = E-TS- ]T ViN*, (2) 

i=u,d,s 



from which the relevant equations of state for the entropy S, the pressure P 
and the particle number can be calculated as usually (the expressions are 
given in Sec. IV of Ref. [20]). So we take the temperature T, the volume V 
and the chemical potential of the i-quark (/ij) as the full independent state 
variables, and a grand canonical approach is applied to our model of strong 
interacting matter. It can simulate either a region in the interior of a neu- 
tron star, or a hot and dense fireball created in a heavy-ion collision. Since 
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electrons and positrons are not involved in the strong interaction, we impose 
the condition p e = 0. So we naturally get the chemical equilibrium condition 
li-u — A*d = Pa — that is used along the work. This condition is only valid 
around the CEP, where the temperature is expected to be high enough to 
avoid color condensation, or at small chemical potential where this phenom- 
ena is not yet present. The baryon number susceptibility is the response of 
the baryon number density ps(T, pi) to an infinitesimal variation of the quark 
chemical potential pi [2T] : 




Here, the quark density is pi = Ni/V — % / p 2 dp (n,i(pi, T) — n^p^ T)). Other 
relevant observable, in the context of possible signatures for chiral symmetry 
restoration in the hadron-quark transition and for transition from the hadronic 
matter to the QGP [21,22,23], is the specific heat which is defined by 



C 




T 
V 




\{dN t /dT), t 

(dNi/dpi), 



(4) 



where we have transformed the derivative (dS/ dT) ^. using the formula of the 
Jacobian. In fact, we work in the grand canonical ensemble where (T,V,pi) 
are the set of natural independent variables (still holding iVj and V fixed). 



3 The critical end point of QCD 



To study the influence of explicit chiral symmetry breaking on the location 
of the critical points we vary the current quark masses rrii keeping the other 
model parameters. The phase diagram for the SU(3) NJL model is presented in 
Fig. [1] as a function of ps and T, and considering different cases for the current 
quark masses. The main motivation is the discussion of the critical behavior 
of the system, starting with the location of the CEP. Using physical values of 
the quark masses [18,24]: m u = = 5.5 MeV, m s = 140.7 MeV, this point 
is localized at T CEP = 67.7 MeV and p c B EP = 318.4 MeV (p c B EP = 1.68p )- 
We also verified that, contrarily to what happens in the chiral limit for the 
SU(2) sector where the TCP is found, the NJL model in SU(3), also in the 
chiral limit (m u = ma = m s = 0), does not exhibit a TCP 0: chiral symmetry 
is restored via a first order transition for all baryonic chemical potentials 

2 Both situations are in agreement with what is expected: the chiral phase transition 
at the chiral limit is of second order for Nf = 2 and first order for Nf > 3 [12] . 
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Fig. 1. Left panel: the phase diagram in the SU(3) NJL model. The solid line rep- 
resents the first order phase transition. The size of the critical region is also plotted 
for Xb/Xb 6 = 2(3). Right panel: the phase diagram and the "line" of TCPs for 
m u = md = and different values of m s (the dotted lines are just drawn to guide 
the eye). 

and temperatures (see left panel of Fig. [1]). This pattern of chiral symmetry 
restoration remains for m u = m d 



and m s < mf [25]. In our model we 
found mf lt = 18.3 MeV for m u = m d = 0. When m s > mf lt , at fi B = 0, the 
transition is of second order and, as \i& increases, the line of the second order 
phase transition will end in a first order line at the TCP. Several TCPs are 
plotted for different values of m s in the right panel of Fig. [TJ As m s increases, 
the value of T for this "line" of TCPs decreases as \xb increases getting closer 
to the CEP and, when m s = 140.7 MeV, it starts to move away from the 
CEP. The TCP for m s = 140.7 MeV is the closest to the CEP and is located 
at [i T B CP = 265.9 MeV and T TCP = 100.5 MeV. If we choose m u = m d ^ 0, 
instead of second order transition we have a smooth crossover for all the values 
of m« and the "line" of TCPs becomes a "line" of CEPs. 



4 Behaviour of xb and C in the vicinity of the CEP and their 
critical exponents 



A bound to the size of the critical region around the CEP can be found by 
calculating the baryon number susceptibility, the specific heat and their critical 
behaviours. If the critical region of the CEP is small, it is expected that most 
of the fluctuations associated with the CEP will come from the mean field 
region around the critical region |10j . 

In the left panel of Fig. [2] is plotted the baryon number density for three 
different temperatures around the CEP. For temperatures below T CEP we 
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Fig. 2. Baryonic density (left panel) and baryon number susceptibility (right panel) 
as function of \xb for different temperatures around the CEP: T CEP = 67.7 MeV 
and T = T CEP ± 10 MeV. 

have a first order phase transition and, consequently, xb has a discontinuity 
(right panel of Fig. [2]). For T = T CEP the slope of the baryon number density 
tends to infinity at \ib — Hb EP which implies a diverging susceptibility (this 
behaviour was found in [TOfTl] using different models in the SU(2) sector). 
For temperatures above T CEP , in the crossover region, the discontinuity of 
Xb disappears at the transition line, and the density changes gradually in a 
continuous way as we can see in the right panel of Fig. [2l A similar behavior 
is found for the specific heat for three different chemical potentials around the 
CEP, as we can observe from Fig. [3j 

As we have already seen, several thermodynamic quantities diverge at the 
CEP. In order to make this statement more precise, we will focus on the values 
of a set of indices, the so-called critical exponents, which describe the behavior 
near the critical point of various quantities of interest (in our case e and a are 
the critical exponents of xb and C, respectively). The motivation for this study 
arises from fundamental phase transition considerations, and thus transcend 
any particular system. These critical exponents will be determined by finding 
two directions, temperature-like and magnetic-field-like, in the T — /ig plane 
near the CEP, because, as pointed out in [26], the form of the divergence 
depends on the route which is chosen to approach the critical end point. 

Starting with the baryon number susceptibility, if the path chosen is asymptot- 
ically parallel to the first order transition line, the divergence of xb scales with 
an exponent jb- In the mean-field approximation it is expected 7b = 1 for 
this path. For any other path not parallel to the first order line, the divergence 
scales with the exponent e = 1 — 1/5. Once in the mean-field approximation 
5 = 3, we will have e = 2/3 and 7b > e is verified. The last condition is respon- 
sible for the elongation of the critical region, xb being enhanced in the direc- 
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Fig. 3. Specific heat as a function of T for different values of \ib around the CEP 
(pCEP = 318 A Me y and /JB = ^cep ± 1Q Me y) 

tion parallel to the first order transition line (see Fig. [1]). To estimate the crit- 
ical region around the CEP we can calculate the dimensionless ratio Xb/Xb 6 
where Xb^ is obtained taking the chiral limit m u = = m s = 0. Left panel 
of Fig. [H shows a contour plot for two fixed ratios {xb/Xb^ = 2.0(3.0)) in 
the phase diagram around the CEP where we confirm the elongation, in the 
direction parallel to the first-order transition line, of the region where xb is 
enhanced. 



For the baryon number susceptibility we will start with a path parallel to the 
/ie-axis in the (T, /z^-plane, from lower hb towards the critical n E [ EP = 318.4 
MeV, at fixed temperature T CEP = 67.7 MeV. In Fig. @] we plot xb as a 
function of \xb close to the CEP. Using a linear logarithmic fit 

ln%B = -eln|/x B -/4? BP | + c i> ( 5 ) 



where the term c\ is independent of /ib, we obtain e = 0.67 ± 0.01, which is 
consistent with the mean field theory prediction e = 2/3. 

Once there is no reason why the critical exponent should be equal for both 
regions, below and above /^ EP , we also study the baryon number susceptibility 
from higher jj,B towards the critical fi^ EP - The logarithmic fit used now is 
Iio-Xb = — e' In \hb — P \ + c l- O ur result shows that e' = 0.68 ± 0.01 which 
is very near the value of e. This means that the size of the region we observe 
is approximately the same independently of the direction we choose for the 
path parallel to the /xg-axis. 

Paying attention to the specific heat around the CEP, we have used a path 
parallel to the T-axis in the (T, /is)-plane from lower (higher) T towards the 
critical T CEP = 67.7 MeV at fixed yP B EP = 318.4 MeV. In Fig. 1 we plot C as 



7 




Fig. 4. Left panel: baryon number susceptibility as a function of \hb — Ms I a ^ 
fixed temperature T = 67.7 MeV. Right panel: specific heat as a function of 
\T - T CEP \ at fixed chemical potential yf^ EP = 318.4 MeV. 

a function of T close to the CEP in a logarithmic scale. Now, we see that for 
the region T < T we have a nontrivial critical exponent a = 0.61 ±0.010. 
This value of a is not in agreement with what was suggested by universality 
arguments in [TO] : it is expected that xb and C should be essentially the 
same near the TCP and the CEP which implies a = e = 2/3. One possible 
interpretation of this result could be the effect of the hidden TCP on the CEP, 
as it was already seen in Refs. [T0l[l~4"] . and which influence could be stronger 
in the specific heat rather than in the baryon number susceptibility. However, 
this explanation is not valid in the framework of the NJL model once, at the 
TCP, the value of a is already not consistent with the respective mean field 
value. 

To support the last statement, let us analyze the behavior of the specific heat 
around the TCP when m s = 140.7 MeV (T TCP = 100.5 MeV, n T B CP = 265.9), 
the nearest TCP from the CEP (see right panel of Fig. [I]). Using a path 
parallel to the T-axis from lower T towards the critical T TCP , at fixed chemical 
potential /i^ CP , we find a = 0.45 ± 0.01. This result, in spite of being close, 
is not in agreement with the respective mean field value (a = 1/2). We notice 
that the inconsistency with the mean field values only occurs for the specific 
heat. In fact, it is found that the critical exponent for \b (jb once we are in 
the TCP) whose value, 7# = 0.50 ± 0.02, is in agreement with the respective 
mean field value (js = 1/2). 

Nevertheless we observe that the values of a in the TCP and in the CEP are 
consistent within the NJL model. We also stress that the universality argu- 

3 We use the linear logarithmic fit In C = — a In \T — T CEP \ + C2 where the term C2 
is independent of T . 
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ments are so general that they give no quantitative results and, due to the lack 
of information from the lattice simulations, they should be confronted with 
model calculations. The eventual difference between the values of the C and 
Xb critical exponents can be interesting in heavy-ion collisions experiments. 
Finally, for the region T > T CEP the critical exponent is a' = 0.67 ± 0.01 
which is compatible with the value of e'. This means that the specific heat is 
sensitive to the way we approach the CEP. 



5 Summary and conclusions 

We have analyzed the phase diagram in the SU(3) NJL which reproduces the 
essential features of QCD: a first order phase transition for low temperatures 
and the existence of the CEP. Contrarily to what happens in the SU(2) sector, 
in the chiral limit (m u = = m s = 0) we do not find a TCP in the NJL 
model, which agrees with what is expected: the chiral phase transition at = 
is second order for Nf = 2 and first order for Nf > 3. When m u = = 
and m s > m^"* (mf = 18.3 MeV) the transition is second order ending in 
a first order line at the TCP. As m s increases we have a "line" of TCPs. For 
m u = 7^ there is a crossover for all the values of m s and the "line" of 
TCPs becomes a "line" of CEPs. The location of the CEP depends strongly 
of the strange quark mass. Around the CEP we have studied the baryon 
number susceptibility and the specific heat which are related with event-by- 
event fluctuations of /ie or T in heavy- ion collisions. In the NJL model, for 
Xb, we conclude that the obtained critical exponents are consistent with the 
mean field values e = e' = 2/3 (the NJL model only produces mean field 
behaviors). From our study of the critical exponent for the specific heat, we 
conclude that a is different from e. More relevant information about the CEP 
can be obtained from the spectral functions and the isentropic trajectories in 
its vicinity. This work is in progress. 
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